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Q\ . Abstract 

We argue that the spectrum of the QCD Dirac operator near zero virtuality can be 
^' described by random matrix theory. As in the case of classical random matrix ensem- 



bles of Dyson we have three different cases: the chiral orthogonal ensemble (chGOE), 
the chiral unitary ensemble (chGUE) and the chiral symplectic ensemble (chGSE). They 
correspond to gauge groups SU{2) in the fundamental representation, SU{Nc), Ac > 3 in 
the fundamental representation, and gauge groups for all Ac in the adjoint representation, 
respectively. The joint probability density reproduces Leutwyler-Smilga sum rules. 
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According to the Banks-Casher formula |jl|, the spectrum of the Dirac operator near zero 
virtuahty is directly connected with a nonzero value of the chiral condensate, the order 
parameter for the chiral phase transition. This suggests that the spectrum in this region 
plays an important role in understanding the mechanism of chiral symmetry breaking. 
Recently, it was shown that a nonzero value of the chiral condensate leads to the existence 
of sum rules for inverse powers of the eigenvalues of the Dirac operator p[. These sum 
rules are only sensitive to the spectrum near zero virtuality and can be expressed into 
microscopic spectral correlation functions, which measure correlations on the order of a 
finite number of average level spacings [0, Q • As is well known from the study of chaotic 
systems p, ^, such correlations are independent of the details of the interactions and 
can be described by random matrix theory. Indeed, we |^, ^ have shown that for of 
a complex Dirac operator all Leutwyler-Smilga sum rules follow from a chiral random 
matrix theory. This led to the claim that the microscopic spectral correlation functions 
of the Dirac operator are universal. 

More than three decades ago Dyson [^ found three distinct types of random matrix en- 
sembles: the gaussian orthogonal ensemble (GOE), the gaussian unitary ensemble (GUE) 
and the gaussian symplectic ensemble (GSE), corresponding to real, complex and quater- 
nion matrix elements. The random matrix theory discussed in [Q, ^] has complex matrix 
elements. For that reason it has been named it the chiral unitary ensemble (chGUE). 
This raises the question of what are the chiral analogues of the GOE and the GSE, and 
what is the structure of the corresponding Dirac operator. The answer to this question 
will be given in this letter. We also will derive the joint eigenvalue density of the random 
matrix ensembles and present the result for the simplest Leutwyler-Smilga sum rule. 

The Euclidean Dirac operator in QCD is defined by 

D = i-fd + -fA, (1) 

where A is an SU{Nc) valued gauge field {Nc is the number of colors). Because this 
operator anti-commutes with 75, in a chiral basis it reduces to the following block structure 

In general (for Nc > 3) , ii A is in the the fundamental representation, the matrix elements 
of the Dirac operator, i.e. Tij, are complex. This defines the first family of Dirac operators. 



However, in the case of SU{2) we have an additional symmetry p[, which is specific 
to this group: 

[C-Sir,D] = 0, (3) 

where C is the charge conjugation operator (7* = — C7^C^^), and K denotes the complex 
conjugation operator. This symmetry operator has the property that 

{C-^T^Kf = 1. (4) 

As is well known from the analysis of the time-reversal operator in random matrix theory 
0, this property allows us to choose a basis in which the matrix elements of the Dirac 
operator in (2) are real and T^ = T. This provides us with the second family of Dirac 
operators. 

The third family of gauge theories are those with the fermions in the adjoint repre- 
sentation. The Dirac operator is given by 

Dm = iiddui + f"^lAm, (5) 

where f^^"^ are the structure constants of the gauge group. As was noted in 0, in this 
case the covariant derivative is real, and the Dirac operator is invariant under charge 
conjugation 

[D,C-^K] = Q. (6) 

Because C*C^^ = — 1, one can easily derive that 

[C-'K)' = -l. (7) 



This implies that each eigenvalue of the Dirac operator is doubly degenerate with linearly 
independent eigenfunctions p| given by 

(t)x and C~^K(t)x. (8) 

As follows from a discussion by Dyson [^], in this case the Dirac operator can be diago- 
nalized by a symplectic transformation. Or, in other words, the matrix elements of T can 
be regrouped into real quaternions, and T}j = Tji (quaternion conjugation is denoted by 
a bar). 



The QCD partition function for Nf flavors with masses rrif {rrif -^ 0) in the sector 
with u zero modes is defined by 

Nf 

Z?'''' = {Il n(^n + ^?K)mA), (9) 

/=1A„>0 

where the average < ■ ■ ■ >s^{A) is over gauge field configurations with u fermionic zero 
modes weighted by the gauge field action S^{A). The product is over all eigenvalues of 
the Dirac operator. For fermions in the adjoint representation, the doubly degenerate 
eigenvalues count only once (Majorana fermions) 0. Relevant observables are obtained 
by differentiation with respect to the masses. The distribution of the eigenvalues of the 
Dirac operator is induced by the fluctuations of the gauge field. It is our claim that the 
correlations between eigenvalues near zero virtuality, i.e. only a finite number of level 
spacings away from zero, do not depend on details of the interaction. 

The basic underlying idea of random matrix theory is that correlations between eigen- 
values on the scale of one eigenvalue are only determined by the symmetries of the system 
and do not depend on the detailed dynamics. Therefore, the spectral density measured 
in units of the average spectral density near zero virtuality is a universal quantity that 
can be described by a random matrix theory that reflects only on the symmetries of the 
Dirac operator. The relevant random matrix theory in the sector with u zero modes is 

Z,,. = J VTP,{T) n det ( "^^^ '^^ ) , (10) 

where T has the symmetries of the corresponding Dirac operator and the masses are in 
the chiral limit {m,f -^ 0). For SU{2) in the fundamental representation T is real {(3 = 1), 
for SU{Nc), Nc > 3, the matrix T is complex {/3 = 2), and for fermions in the adjoint 
representation, T is quaternion real (/? = 4). In the latter case the square root of the 
fermion determinant appears in (10). The matrix T is a rectangular n x m, matrix with 
\n — m,\ = V. It can be shown that the matrix in the fermion determinant in (10) has 
exactly v zero eigenvalues (and Z ~ 11/ ^^/)- The distribution function of the matrix 



elements P{T) that is consistent with no additional information input is gaussian [^. In 
a standard normalization we choose 

pm=^w{-^,t>i), (11) 



where the sum is over the non-zero eigenvalues of T. With this choice the average spectral 
density does not depend on P and is given by p(0) = l/vra. This allows us to identify 
a = l/7rp(0) = 1/S, where S denotes the chiral condensate. It should be noted that as 
in the QCD partition function the chiral symmetry is broken isotropically in flavor space. 
In order to derive the joint eigenvalue density we use the eigenvalues and eigenangles of 
T as new integration variables. In each of the three cases, up to a constant, the Jacobian 
of this transformation is given by 

J = U\^l-^yilK'^'-'- (12) 

The derivation of this result will be given elsewhere. At the moment we only remark that 
the total powers of A^ can be obtained on dimensional grounds only. The joint eigenvalue 
density is therefore given by 

k,l k 

The normalization constant is denoted by Cn,i3- In the case of /3 = 4, each of the doubly 
degenerate eigenvalues is counted only once. This is consistent with the fact that fermions 
in the adjoint representation are Majorana fermions (see for a detailed discussion of 
this point). The simplest Leutwyler-Smilga sum rule in the sector with v zero modes can 
be evaluated with the help of Selberg's integral [^ (see [0 for a discussion). The result 
is 

\^AS^/,/8(f + f + ^/-l)' ^''^ 

where the average is with respect to the spectral density p^(Ai, ■ ■ ■ , A„). The total number 
of modes is denoted by A^ = -m + n. This constitutes the final result of this letter. It 
agrees with sum rules obtained by Leutwyler and Smilga for (3 = 2 and /? = 4. This 
sum-rule can also be expressed in terms of the microscopic spectral density defined by 

where p(A) is obtained from (13) by integrating over all eigenvalues except one. The 
microscopic spectral density has been derived for (3 = 1 [l^ and /5 = 2 |Q, and in both 



cases it agrees with numerical results from simulations of gauge field configurations by a 



liquid of instantons [15 



In PI sum rules were derived from the static limit an effective field theory. The above 
discussed triality implies that we have three structurally different effective field theories. 
Two of them were analyzed in B, and the theory for /3 = 1 which involves both baryons 



and mesons, was discussed in |T^. It would be instructive to derive (14) for /3 = 1 in the 
framework of this model. 

In conclusion, we have argued that, depending on the gauge group, the QCD Dirac op- 
erator falls into three different families: SU{2) in the fundamental representation, SU{Nc), 
A'"c > 3 in the fundamental representation and SU{Nc) in the adjoint representation. This 
triality corresponds to real, complex and quaternion matrix elements. Its spectrum near 
zero virtuality reflects only the symmetries of the system and can be described in terms 
of chiral random matrix theory: the chGOE, the chGUE and the chGSE, respectively. 
Sum rules obtained general arguments based on effective field theory |^ are reproduced. 
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